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We calculate the cross sections and final state distributions for the processes e+e^ — >■ 
T{\S){'k'^'k~ , K~ ,ri-K^) near the T(5S') resonance based on the tetraquark hypothesis. This 
framework is used to analyse the data on the T(lS')7r'^7r~ and 'T{1S)K'^ K~ final states [K. F. Chen 
et al. (Belle Cohaboration), Phys. Rev. Lett. 100, 112001 (2008); I. Adachi et al. (Belle Collabo- 
ration), ,arXiv:0808. 2445 , yielding good fits. Dimeson invariant mass spectra in these processes are 
shown to be dominated by the corresponding light scalar and tensor states. The resulting correla- 
tions among the cross sections are worked out. We also predict a{e^e~ T (1S)K^ K~) / a^e^ e~ — >■ 
T(IS') A'^A''') — 1/4. These features provide crucial tests of the tetraquark framework and can be 
searched for in the currently available and forthcoming data from the B factories. 



PACS numbers: 14.40.Rt, 14.40.Pq, 13.66.Bc 



The anomalously large production cross sections 
for e+er ^ T {IS)tt+ tt' ,T {2S)n+ tt- and e+e~ 
T(3S')7r+7r^ measured between -^s = 10.83 GeV and 
11.02 GeV by the Belle Collaboration [1 Q at KEK 
do not agree well with the lincshape and production 
rates for the conventional hh state T(10860) (also called 
T(5S')). A fit to the measured production cross sec- 
tions using a Breit-Wigner resonance shape yielded a 
peak mass of [10888.4t2;g(stat) ± 1.2(syst)] MeV and a 
width of [30.7t7;o(stat)±3.1(syst)] McV for the observed 
state, henceforth called Yh(10890) 0. More data arc 
required to understand the resonance structure in this 
region. In Q, a dynamical model was developed to ex- 
plain the Belle data for the final states T(lS')7r+7r^ and 
T(25')7r+7r^ in terms of the production and decays of 
the states ^jt,,;//!], which are linear superposition of the 
J^*^ = 1 hidden hh tetraquark states Yjbu] = 
and F[f,d] = The mass difference, estimated as 

M{Y^b,h\) - M{Y^t,j{) = (5.6 ± 2.8) MeV 0, i, was ig- 
nored and the mass-degenerate states ^jb,;] and i^[b,/i] were 
identified with the Fb(10890). This model described the 
distributions in the dipion invariant mass and the helicity 
angle measured by Belle [ij well and offered an explana- 
tion of the rates in terms of the Zweig-allowed transitions 
y[b.i/h] (T(lS'),T(25'),T(35'))7r+7r-. In the case of 
A (3872), the phase space for the decay A (3872) DD* 
is highly constrained, yielding a very small binding en- 
ergy: So ~ Mx - Md - Md- = -0.25 ± 0.40 McV. 
This led to suggestions that A (3872) could be a hadron 
molecule. In the case of yb(10890), the phase space is 
large enough to allow the decays r;,(10890) B^*^B^*\ 
and consequently the decay width is large, Vy^ — 30 
MeV, which is of the same order as the total width of 
the T(5S'). Hence, (10890), in all likehhood, is not a 
hadron molecule. While credible, the interpretation of 
yb(10890) in terms of the hb tetraquark states requires 
further experimental scrutiny. It is the aim of this Letter 
to provide some definitive tests to confirm or rule out the 



tetraquark interpretation of the Belle data. 

To that end, we further develop the tetraquark formal- 
ism for the processes 

e+ + e- ^ Y^{q) ^ T(l^)(p) -f- P{ki) + P'ik^) , (1) 

where PP' stands for the pseudoscalar-meson pairs 
TT+TT^, A'+AT" and rjir'^, and q, p, ki and k2 are the 
momenta of Yt, T(IS'); P and P', respectively. We ne- 
glect other background processes, based on prior data on 
the dipionic transitions involving higher T{nS) to lower 
T{mS) (to < n) states Following Q, the low mass 
scalar 0++ hadrons a or /o(600), /o(980) and a[J(980) 
(the upper index indicates the /s = component of the 
iso-triplet ao), which enter as intermediate states in the 
processes n ^ T(IS') + [a, /o(980), • • •] ^ T[1S)PP', 
are assumed to be tetraquark states. These intermedi- 
ate = 0++ states together with the J^^ = 2++ 
state /2(1270) provide the dominant resonating part of 
the amplitudes for the processes considered in this work. 
These resonances are labeled cr, /o, ftg and /2 hence- 
forth. We determine the coupling constants involving 
these light tetraquark states and the mesons PP' from 
the known decays from the PDG 0] and data from the 
E791 [3, the BES @, the Crystal Barrel (CB) and 



the KLOE Collaborations [UJ, II4I , adopting the Flatte 



model [ist for the cr, /o and Cq couplings to take into ac- 
count threshold effects. The non-resonating continuum 
contributions are parameterized in terms of two a pri- 
ori unknown constants [14| . With this formalism, we 



analyze the invariant-mass Mppi and the cos 9 spectra, 
where Mppi ~ {ki + /s2)^ and is the angle between the 
momenta of Yf, and P in the PP' rest frame. 

The theoretical framework described here provides 
good fits of the Belle data on the invariant dipion mass 
spectrum and cosO distribution in the process e^e~ — )■ 
T(lS')7r+7r" and the ratio crT(is)if+if-/CTT(is)7r+7r- , 
with iyT(is)PP' being the cross section cr(e"'"e^ — > 
T(IS')PP'). We present the invariant mass distribu- 
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tions for the K'^K~ and -q-K^ mesons in the processes 
e^e~ — >■ T(lS')(if+i^^, yyTT*^), which are dominated by 
the respective J^'^ = 0++ resonances. The result- 
ing correlations among crx(i5)7r+7r- , o-t{is)k+k- and 
CTy(i5)^^o are worked out. Constraining these corre- 
lations from the existing data on the first two pro- 
cesses, we predict crT(is)r)iro/f T(is)ir+7r- ■ We also pre- 
dict tTx(is)if+/f-/o'T(is)KO/?o = 1/4, reflecting the ratio 
Qfbu]/Qfbd] '^itli Qlbn] = 1/3 and = -2/3 being the 
effective electric charges for the constituent diquarks of 
Y[bu] and Y[m] , respectively. 

We start by defining the tetraquark states in the 
isospin basis, with the two isospin components = 
(Y[bu] + Y[M])/V2 and Yf^ = {Y[i,u] - Y[bd])/V^ for isospin 
/ = and / = 1, respectively. The two mass eigenstates 
Y[bfi and Y\f,.h] are identified with Y[h„,i and Yf^rfi, as the 



f[b,h] 



[bu] 



'[bd], 



mixings between them is small. We ignore the mass dif- 
ference and also the isospin breaking effects except for 
the production processes e+e^ — > Y"/ hereafter. 

We calculate the decay amplitude as a sum of the Breit- 
Wigner resonances and non-resonating continuum contri- 
butions, with the latter adopted from IJ]. The differen- 
tial cross section is then written as 



(IS)PP' 



dMpp' dcos0 ~ 3847r3s Mpp> [(s - m2,J2 4. m|.^r2,J 
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where A(a;, y,z) = {x — y — z)^ — 4yz, g° and p° are the 
energies of the Yb and T(15') in the PP' rest frame, re- 
spectively, Ty^ is the decay width of Yb, and my^, rnr, 
mp and mp' are the masses of Yf,, T(IS'), P and P', re- 
spectively. We take my, = 10.89 GeV and Ty^ = 30 MeV 
in the numerical analysis below. A detailed derivation of 
the above formula will be presented in [l5| . 

Each PP' channel receives specified contributions de- 
pending on the isospin of PP' and the kinematically 
allowed region for the invariant mass Mppi e [mp + 
mpi ^ y/s — mx(is)]. The S-wave amplitude for the PP' 
system, S, and the D-wave amplitudes, P, V and P", 
are the sums over possible isospin states 



M=^Mi for M=S,V, V, V", 



(3) 



where / = for tt+tt", / = 0, 1 for K+K', and / = 1 
for TyTr", since the T(15) is an isospin state, and the 



following resonances contribute to each process: 
a, fo and /2 for tt+tt". 



/o 



.0 



and /2 for K+K~ 
for TyTT^ . 



(4) 



The 7 = amplitudes are given by the combinations of 
the resonance amplitudes, and TMq^, and the non- 
rcsonating continuum amplitudes, and M.^'-^: 



Vo^\k\^Mt 



v. 



2?o , 



2gV 



0^0 



(5) 



(9 -P) 

where S runs over possible / = scalar resonances in 
Eq. (HI, and |fe| is the magnitude of the three momen- 
tum of P^'^ in the PP' rest frame. Similarly, the / = 1 
amplitudes are given by 



5l : 



5e+e 
^0, 



(fci • fc2)A^?' 
5e+ 



V'[ ^ f^^Al^c^, (6) 



where the dimensionless couplings g^+e-Y^' are defined 
through the Lagrangian C = E/=o,i 9e+e-Y^'Ybt, (e7^e), 
and the ratio is given by Se+e-y.i/ffe+e-y," = (Q[H - 

Q[bd])/iQ[bu] + Q[bd]) = -3. 

To calculate the production cross sections, we derive 
the corresponding Van Royen-Weisskopf formula for the 
Icptonic decay widths of the tetraquark states made up 
of point-like diquarks: 



^{Y[bu/bd] e+e ) = 



2Aa'^\Qlbu/bd]\ 



i?lV(0) , (7) 



where a is the fine-structure constant, the parameter 
K takes into account differing sizes of the tetraquarks 
compared to the standard bottomonia as well as the 
hadronic size of the diquarks, with k < 1 anticipated, 
and |i?(V(0)|2 = 2.067 GeV^ ^ is the square of the 
derivative of the radial wave function for Xbi^P) taken at 
the origin. Hence, the leptonic widths of the tetraquark 
states are estimated as 

r(Y[M] ^ e+e-) = ATiY[bu] ^ e+e") « 83 eV , (8) 

which are substantially smaller than the Icptonic width of 
the T(5S') 0. Combining the knowledge of the produc- 
tion process and the measured cross section for e+e" — 
T(1S')tt+tt", we estimate T{Y[bu/bd] ^ T(1S')tt+tt") to 
be 0(1) MeV. 

The continuum amplitudes in Eq. ( El) a re written in 
terms of the two form factors A and B [ij, [l5| as 



2A' 
B' 



(fci • ^2) + 



B' 3(g")'fc;fcg- |gp|fc|' 
fpfpi 3s 



fph 



(9) 
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where the primed quantities here, and later, are rescaled 
as A' = Ag^+^-yo and B' = Bg^+^-yo, fpio is the decay 
constant of P''\ and \q\, fc° and are the magnitude of 
the three momentum of Yt and the energies of P and P' 
in the PP' rest frame, respectively. Using SU(3) symme- 
try for the form factors A and B in Eq. ©, but not for 
the pseudoscalar meson masses and coupling constants, 
we assume the relations M^'"'^'" {T{1S)K'^K^) = 
(v^/2)7Wj^''^(T(15)7r+^-), Ml^'^^{T{lS)K+K-) ^ 
Ml'^'^^iT{lS)K+K-) and 7W}^'^^(T(15)?77rO) = 
V2 MI^'"^^ {TilS)K+ K-). 

The resonant contributions are expressed by the Brcit- 
Wigner formula: 



i iriRTR 



(10) 



where / = for i? = cr, /o and /2, and 7=1 for 
R = Uq. The couplings for the scalar resonances S arc de- 
fined through the Lagrangian£ = gspp'{dfiP){d^^P') S+ 
9YbT{is)sYbfiT'^S, while those for the /2 are defined via 
C = 2gf,pp,id^P){d,P')fr +9nTiis)f,Yb,,T,fr ■ The 
couplings gRppi and gyif(is)R. have mass dimensions —1 
and 1, respectively. For the ct, /o and Og, we adopt the 
Flattc model ^ 

'm-a^a = faTTTrPTTTr , ^fo^fo — //oTTTrPTTTr + ffgKRPKK i 
rriaO^a'^ = fa«n^Pr,n + IIohRPkR (H) 

with the phase space factor pab = [(1 — ("nia — 
rubf /Mlp,){l-{ma+mbf /Mlp,)Y/'^ , where the Flatte 
couplings fspp' are related to the vertex couplings gspp' 
entering in Eq. (jlOp via 



gspp'ih ■ k2) = 4,\^ fspp'- (12) 
The couplings in Eq. (jl2p arc defined for an exclusive 



final state. Summing over final states, e.g., ff^T^T, = 
//07I-+7I- +//o7r07rO' obtain the isospin relations /sxtt = 

fsn + Tv- 7 fsKK = V^fsK+K- and fsr,TT = /sj)7rO- 

For the a meson, we extract the coupling 17^^+^- from 
the E791 data Q: = 26.7 GeV~^ with = 478 

MeV, yielding the Flatte coupling f^,,.^ = 437 McV. For 
the /o and Oq mesons, wc adopt the masses and the Flattc 
couplings measured by the BES 0] and CB [l^l Collabo- 
rations (the corresponding couplings from KLOE jlllll2| 
arc shown in the parentheses): 

=965(984), =406(349), 4^,^,^ = 833(869), 
m,o =982(983), /„o^, = 324(398), /.o,,.,^ = 329(429) (13) 

in units of McV. Furthermore, we extract the couplings 
for the /2 meson from r(/2 PP') = g'j^pp,m^j^{l — 
4m|,/TO^^)^/^/(4807r) for nip = mpi, where the data for 
r(/2 ^ ^tt) = (3/2) r(/2 ^ n+TT-) and r(/2 KK) = 
2r(/2 K+K~), and m/^ = 1275 McV arc taken from 
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FIG. 1. Fit results (a) of the M^+^- distribution and (b) of 
the cosS distribution for e^e^ — > Y;, T(lS')7r"'"7r^, normal- 
ized by the measured cross section. The histograms (green 
bars) represent the fit results, while the crosses are the Belle 
data [31 . The resonant contributions from the a and /o(980) 
(left red curve) and the /2(1270) (right black curve) are also 
indicated in (a). 



TABLE I. Best fit parameters, yielding x^/d.o.f. = 21.5/15, 
where A' and B' are dimensionless, ffyo-i-(is)/,, ^''id (;yo-i-(is)/2 
are given in units of MeV, and the angles are in units of rad. 
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0.000079 


-0.00020 0.318 0.439 


0.36 


-2.76 


-0.46 



PDG |6| . The other inputs for the pseudo-scalar mesons 
and the T(IS') are also taken from PDG. 

Having detailed our dynamical model, we now per- 
form a simultaneous fit to the binned T(lS')7r+7r^ 
data for the Af^+^- and cosd distributions measured 
by Belle at ^/s = 10.87 GeV [l[, normalizing them 
by the measured cross section: da^+^- /dM.^+^- and 
da„+„-/dcos0, where = crr{is)yr+yr- /<7rfis}7T+7T- 

with erf = 1-61 ± 0.16 pb Q. With SU(3) sym- 

metry for the Y^T{1S)R couplings entering in Eq. ()10p . 
i.e., setting gyar(is}a = 5y,"t(is)/o. ^e have 7 free pa- 
rameters: 

A', B' , 5yOx(iS)/o> 9Y°T{lS)f2^ "^/f" ^S^' ^1^) 

Hence the number of degrees of freedom (d.o.f.) in the 
fit is (12 -I- 10) — 7 = 15. To make predictions for the 
T{1S)K+K- and T{lS)r]n° modes, we assume SU(3) 
for the couplings and the phases, i.e., gY°T(is)fo — 
9Yi^r{is)a° and (pf^ = ip^^o. We also assume that there 
is no phase difference between the two continuum ampli- 
tudes MI'^'^^ and MI^'^^. 

With these inputs, we have performed a large num- 
ber of fits (typically 0(5000)) of the Belle data with the 
tetraquark theory predictions. The resultant best fit is 
fairly good, with xVd.o.f. = 21.5/15 for the BES and CB 
input in Eq. (fT3|) . which corresponds to a p- value of 0.12. 
The corresponding best fit using the KLOE data is very 
similar, having a x^/d.o.f. = 21.9/15, yielding a p-value 
of 0.11. The best fit using the BES and CB data is pre- 
sented in Fig. [T] and the corresponding fit values of the 
parameters arc listed in Table HI Further details about 



4 




Mk+k- [GeV] M^7,o [GeV] Sa'+A'- 

FIG. 2. Predictions (a) of the M^+j^- distribution for e'^e~ — > Vt, — >• T(1S)K^ K~ , (b) of the A/^^o distribution for e"'"e~ — >■ 
Yb T(lS)r;7r'' and (c) of the correlation between the cross sections of T{1S)K'^ K~ and T(lS')777r'', normahzed by the 
measured cross section for the T(lS')7r"^7r~ mode. In (a) and (b), the dotted (solid) curves show the dimeson invariant mass 
spectra from the resonant (total) contribution. In (c), the red dots represent predictions from our fit solutions satisfying 
X /d.o.f. < 1.6. The shaded (green) band shows the current Belle measurement 5^+^^- 

= 0.11^^3 li- 



the correlations among the parameters and the cross sec- 
tions will be presented in a forthcoming paper [l5| . Con- 
cluding the discussion of the final state T(15')7r+7r^, we 
note that the resonance contribution represented by the 
left red curve (S-wave from a and /o) and the right black 
curve (D-wave from /2) in Fig.[T](a) dominate the M^^+j^- 
spectrum, supporting our dynamical model in the decay 
Yb — >■ T(lS')7r+7r~. Sufficient data may provide enough 
statistics to undertake an analysis in the end-region of 
A/jr+TT- to probe the angular distribution of /2 — > tt+tt^. 

The normalized Mj^+j^- and M^j^^o distributions, cal- 
culated with the best-fit parameters in TableHl are shown 
in Fig. [2] (a) and Fig. [2] (b), respectively. In these fig- 
ures, the dotted (solid) curves show the dimeson invari- 
ant mass spectra from the resonant (total) contribution. 
Since these spectra are dominated by the scalars /o + Qq 
and Oq, respectively, there is a strong correlation between 
the two cross sections. This is shown in Fig. [2](c), where 
we have plotted the normalized cross sections ct^+k- 
and cr^Tr" resulting from our fits (dotted points) which 
all satisfy x^/d.o.f. < 1.6. The current Belle measure- 
ment CT/f+if- — 0-llio 03 Q shown as a shaded (green) 
band on this figure. Our model is in agreement with the 
Belle measurement, though there is a tendency in the 
fits to yield larger value for ^k+k- ■ Our predictions will 
be further tested as and when the cross section ct,,^o is 
measured. Noticing that we have neglected the SU(3)- 
breaking effects, we predict 1.0 < ct^ttO ^ 2.0. 

Finally, we note that the states T{1S)K^ K~ and 
T{\S)K'^K^ are produced by the underlying mechanism 
e+e- ^ Y^bu] -> T{1S)K+K- and e+e" ^ ^[^1 ^ 
T{1S)K'^K^ . Hence, a firm prediction is 

'^r{lS)K+K- _ Qffen] _ 1 ^^^^ 
^r(lS)KOK° Qfbd] 4 

This relation is valid under the assumption that the di- 
quarks are point-like. In terms of the mass eigenstates, 
we predict (Tr{is)K+K- = crT{is)KsKs- 



The distributions, cross sections, and correlations pre- 
sented here arc crucial tests of the underlying tetraquark 
hypothesis in the bb sector and go well beyond what has 
been proposed in the literature to understand the nature 
of the yf,(10890) state 0]. They will be scrutinized soon 
in the existing and the forthcoming Belle data to which 
we look forward. 
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